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A WASSERSTEIN GRADIENT FLOW APPROACH TO 
POISSON-NERNST-PLANCK EQUATIONS 

DAVID KINDERLEHRER, LEONARD MONSAINGEON, AND XIANG XU 


Abstract. The Poisson-Nernst-Planck system of equations used to model ionic 
transport is interpreted as a gradient flow for the Wasserstein distance and a free 
energy in the space of probability measures with finite second moment. A varia¬ 
tional scheme is then set up and is the starting point of the construction of global 
weak solutions in a unified framework for the cases of both linear and nonlinear 
diffusion. The proof of the main results relies on the derivation of additional es¬ 
timates based on the flow interchange technique developed by Matthes, McCann, 
and Savare in m- 
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1. Introduction 

The Poisson-Nernst-Planck (PNP) system of equations |18[ is the principal 
description of ionic transport of several interacting species. It has been applied in 
a number of contexts ranging from electrical storage devices to molecular biology, 
at times coupled to Navier-Stokes or other systems. The basic system is to find 
ri(f, x) > 0 and v{t, x) > 0 satisfying 

dtu = Au"^ + div (uV {U + Ip)) , 

dtv = Av'^ + div {vV {V - Ip)), t > 0, x d > 3, (1.1) 

—A^p = u — V, 

for some suitable initial conditions u\t=o = and v\t=o = The unknowns u and 
V represent the density of some positively and negatively charged particles. Here 
m > 1 is a chosen fixed nonlinear diffusion exponent. Note that dEU) formally 
preserves the mass 

/ u{t,x)dx= [ u°(x)dx and / v{t,x)dx= [ v°{x)dx 

jRd- jR<i Js.d- jR<i 

for all t > 0, which physically represents the conservation of total charge of the 
individual species. For initial masses u^{x)dx = v^(x)dx, a simple rescaling of 
time and space allows normalization of masses to unity, and without loss of generality 
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consider u{t,x),v{t,x) to be probability densities. The external potentials U{x) 
and V{x) are prescribed and sufficiently smooth. 'ip{x) from the Gauss Law is the 
self-consistent electrostatic potential created by the two charge carriers according to 
the last equation in dEB. The first two equations in (11.11) are called Nernst-Planck 
equations and describe electro-diffusion and electrophoresis according to the Tick 
and Kohlrausch laws, respectively, while the last equation in (jl.ip corresponds to 
the electrostatic Poisson law. The boundary condition for this coupling equation 
will always be understood in the sense of the Newtonian potential: we shall always 
implicitly write 

— Aij; = u — V <(=> ijj = G * {u — v) = {—A)~^{u — v), (1-2) 


where 


G{x) := 


1 


\2-d 


d{d — 2)uid 


is the fundamental solution of —A in and Ud the volume of the unit ball in 
(for d >3). 

fn some PNP models, an extra background doping profile G{x) is considered 
resulting in the modified potential —Aip' = u — v + G. With suitable assumptions 
on G{x) this can be easily eliminated replacing the external potentials U hy U + il^c 
and V hy V — where ■i/’c = (— 

There is a vast literature on well-posedness and long time behavior of the system 
dlT]). We refer to m n [Ml [n] and references therein for bounded domains, and 
[aiiiiTi [Ml EH for the whole space problem. Different from these papers, our 
contribution is to show that the system of equations governing drift, diffusion 

and reaction of charged species, possesses a gradient flow structure as viewed on the 
metric space of probability measures on endowed with the quadratic Wasserstein 
distance, in essence the weak-* topology (see also |28 1 131 1 136] for various discussions). 
Therefore, variational methods may be introduced to prove global existence of weak 
solutions (see definition below). Motivation for this in terms of energy dissipation 
is offered below. Also, we provide a unihed framework both for linear m = 1 and 
nonlinear m > 1 diffusions. To the best of our knowledge well-posedness in the 
whole space for m > 1 usually requires either high integrability of the initial data 
or initial gradient regularity. We would like to stress that we need here no such 
hypotheses and that our result merely requires some low initial integrability, defined 
in terms of the diffusion exponent m > 1 and dimension d> 3 only (which we think 
is sharp, see discussion after the proof of Proposition 14.41 page I19p . Also, we do not 
need compatibility conditions between m > 1 and d > 3. This is in contrast with 
the Patlak-Keller-Segel chemotaxis models mm. for which critical mass phenomena 
may occur depending on whether m is larger or smaller than the scale-invariant 
exponent m^{d) = 2 — see e.g. |10| for the critical parabolic-parabolic case. This 
is mainly due to the fact that the self-induced drifts are repulsive here, while they are 
self-attractive in the Keller-Segel models, thus leading to aggregation and blow-up 
in finite time. 

It was shown in |16[ (30] that certain scalar diffusion equations can be interpreted 
as gradient flows in metric spaces and the literature concerning this issue is steadily 
growing (see [2] and references therein). It is thus a question of great interest to 
apply such ideas to study systems of equations. In contrast, there are only a few 
examples for systems. For related studies, we refer to |ll[9l[ini[l5lEQl[22l[26l[37|, 
where the energy functional is involved with the Wasserstein metric and existence 



theorems using a minimizing movement scheme for corresponding evolution problemi 
are presented. 

In |16| . the linear Fokker-Planck equation 

dtp = (tA/ 9 + div(/oV<^) 


is regarded as the gradient flow of a free energy consisting of the Boltzmann entropy 
with a potential (p, 


^{P)= i^P +crplogp)dx, 
jR'i 


with respect to the quadratic Wasserstein metric. There may be many Lyapunov 
functions associated to a differential equation. The result of m means that dis¬ 
sipation for the free energy determines the Fokker-Planck Equation. The same 
idea was later employed in |30) and, in addition, to derive long-time asymptotics for 
the Porous Medium Equation (PME). Since the system (II.ip can be viewed as two 
Fokker-Planck equations (when m = 1) or Porous Medium equations (when m > 1) 
in u, V coupled by means of a Poisson kernel, we are motivated to extend these ideas 
to study our coupled system. Inspired by |30| and |35) . we shall in fact discover that 
the PNP system can be seen as a gradient flow driven by the free energy 

£{u,v) := [ (ulogu + vlogv + uU + vV + 


£{u,v) := j (- 


+ 


-|- uU + vV + 


2 

IVV’P 


dx if m = 1, 


(1.3) 


dx if m > 1. 


.m — 1 m — 1 2 

We further motivate this approach informally by discussing the relationship be¬ 
tween the dissipation relation and the weak-* topology in terms of the Wasserstein- 
Rubinstein-Kantorovich distance, or simply the Wasserstein distance. We follow |36) 
and consider for illustration the case m = 1. Set 

(p{u,v) = ulogu + vlogv + uU -\- vV -|- -|V'i/’P, so that 


£{u,v) = / ip{u,v)da 
jRd 


Given a process or an evolution {u{t),v{t)), during an interval (T, T + h) the change 
in energy is 

£{u,v) -I 

T+h Jt 

This, (|1.4I) . is the dissipation equality or inequality and the density of the middle 
term 

-(p{u, u)dx 


—ip(u, v)dxdt = £{u, v) 
dt 


V = - 


dt ‘ 


(1.4) 


(1.5) 


is the dissipation density along the trajectory. Writing 

d , . du du 

we must ascribe a meaning to 

du du 
dt ’ dt 

to render the system dissipative, that is, so that m is positive, 
calculate the terms in (II.6p . Keeping in mind (II.2p . one checks that 

6u' 


( 1 . 6 ) 


To begin we 


= V' and = -V', 








^hich leads to 

ipu{u,v) = logit + U + "0 + 1 and (/?„ = logi; + V — ijj + 1. 

Let us now employ the Poisson-Nernst-Planck equations (HH). Substituting into 
m and integrating by parts gives 


V = 


VU 

-h V(C/ + 'ip) 

u 


u + 


—+ V(F-V') 

V 


u > dx 


Introduce 


w = 


= -( 


Vu 


\ u 


+ V{U + 'il^)] and w = 


= -( 


Vu 


\ V 


+ V{V- 


so that 


ut + div(it;tt) = 0 and vt + div(a;u) = 0. 
We have that 


(1.7) 


rT+h ri+n i- i-i+n p 

/ Vdt = / / \w\^udxdt+ / / \u\‘^vdxdt 

Jt Jt Jt JRd 

where the pairs {u,w), {v,u}) satisfy the continuity equations (ll.7p . This represents 
a trial in dw, the quadratic Wasserstein metric, using the Benamou-Brenier formu¬ 
lation [5], where 

1 pT+h p 


pT+h 


T+h 


dw{u\ 


T+h 


,U\ 


= inf 


liul itdxdt and 


It 

pT+h 


—dw{v\rr,,,,v\rj)‘^=ini l l \uj^vdxdt 

n ^ Jrp Jgd 

taken over all pairs {u,w) and {v,u}) satisfying the continuity equations 
Ut + div(iun) = 0 and vt + div(ti;u) = 0, and 
the initial and terminal conditions. 


The calculation shows that the dissipation relation (11.71) for £ and the PNP system 
are closely related to the Wasserstein distance. We could write, in fact, 

suggesting an implicit scheme which leads to a gradient flow. This is nearly correct. 
As is well known, a factor of 1/2 must be inserted; see below (II.8p . 

We now turn to the precise formulation. Denoting the set of Borel proba¬ 

bility measures on with finite second moments and dw the quadratic Wasserstein 
distance as before, the underlying space will be here {u,v) € P(]R‘^) x and 

will inherit a natural differential structure from that of {V{M.'^),dw) - see section [2] 
for details. The total free energy m is a combination of the well-known internal 
(diffusive entropy) and potential energies for each species, and, although unclear 
at this stage, the coupling Dirichlet energy ^ jVV'pdx falls into the category of 
so-called interaction energies. See |35) for an introduction. 

Following |16| . we shall construct weak solutions 2 : = {u,v) as follows. Given 
suitable initial data = (u^,v^) and some small time step h € (0,1) we first 
construct a discrete sequence solution to the Jordan-Kinderlehrer-Otto or 

JKO implicit scheme 


do) _ m 


= 2 




€ Argmin 
ICxK 




t-.T’)+£(■)} 


( 1 . 8 ) 






Here £{z) = £{u,v) is the total free energy (ll.3p . d? is the (squared) distanci 
on the product space inherited from dw^ and JC C V the set of admissible mini- 
mizers defined later on. As is classical by now, one obtains interpolating solutions 
{zh{t)}h = {uh{t), Zh{t)}h defined for all t>0, piecewise constant in time, and sat¬ 
isfying a coupled system of two Euler-Lagrange equations. We shall then prove that 

as /i —>■ 0 one recovers a weak solution {u{t),v{t)) = z{t) = lim Zhit) of (ll.ip . There 

h^O 

are several challenges in this program. 

In handling the jVV'pdx coupling term, some intrinsic difficulties arise due 
both to the specific Poisson kernel and to the nonscalar setting. First, as neither the 
external potentials nor G are convex the free energy £ is not displacement convex in 
the sense of McCann |27| and we cannot simply apply the standard procedures as 
in [2]. Secondly, due to the singular nature of the kernel G(x) = C both the 
existence of minimizers in (jl.Sp and derivation of the corresponding Euler-Lagrange 
equations become delicate, see in particular the discussions in Proposition 13.11 and 
Proposition l4.4l for details. In order to tackle this issue we used the "flow interchange" 
technique that originates in |25| and was later used in [iniEg to obtain some inte- 
grability improvement and gradient regularity of the minimizers, see Proposition 13.41 
and Proposition 13.51 below. The highlight of the argument is the propagation of ini¬ 
tial LP(]R'^) regularity established in Proposition 13.41 In addition to being technically 
essential here, this propagation of initial regularity allowed us to obtain a natural 
estimate (see Theorems [13]), which to the best of our knowledge was un¬ 
known for the PNP system in the whole space. We also believe that the very same 
argument could be employed for similar problems in order to show propagation of 
initial regularity, which is usually a delicate point in the mass transport framework. 

The rest of the paper is organized as follows. In Section [2] we recall well-known 
facts in optimal transport theory and briefly describe the differential structure of 
the product space. We then formally derive the Wasserstein gradient flow structure 
of the system o and state the main existence results. In Section [3] we study the 
relevant energy functionals, and establish improved regularity of their minimizers. In 
Section m we fix a time step h > 0 small enough and consider the minimizing scheme. 
We obtain approximate discrete solutions {uh,Vh}h and derive the corresponding 
Euler-Lagrange equations. In Section [5] we take the limit h ^ 0 and show that the 

convergence {u,v) = \\m{uh,Vh) is strong enough to retrieve a weak solution. This 
h —>-0 

last section also contains the proof of the main theorems. 

Notation Convention. Unless otherwise specified, (•, •) and • denote inner product 
of elements in V denotes ■p(M'^), and denotes 7^“^(M'^). If clear from the 
context we shall often omit the subscripts m=lorm>l. Ifl<p<oo, we denote 
by p' = the conjugate Lebesgue exponent. 

2. Formal Wasserstein gradient flow 
From now on we assume that the external potentials are quadratic at infinity, i-e 

Cilxp < U{x),V{x) < (72(1 -I- |xp) 

|VU(x)|,|VU(x)| <C3(l + |x|) (2.1) 

ll^b^llL°°(R‘*)i — ^4) 

for some generic positive constants Ci. Note that [/, V need not be uniformly convex 
as is often assumed, so that we allow here multiple wells. Although these assumptions 


&n the potentials could be weakened we assume here strict confinement Ci > 0 for 
the ease of exposition, and we do not seek optimal generality. 

We also introduce the admissible set 


r /Cl := nLilogLi(M'^) ifm = l, 

\ Km ■■= n L^{W^) if m > 1, 


and for reasons that shall become clear later on we shall always consider initial data 


n*’G/C n for some ro > max{m, 2d/(d + 1)}. (2.3) 

Essentially G ensures that the initial energy (ir®, v^) is finite, 

while regularity will ensure that the self-induced drifts riV'I',uV'I' G 

L^{W^) for all times. We believe that r = max{m, 2d/{d + 1)} should be admissible, 
but for technical compactness issues we have to assume here slightly better (M'^) 
integrability for some vq > r arbitrarily close. See the proof of Theorem [¥] later on 
for details. 

For p,u,v > 0 let us define the usual Boltzmann entropy 



plogpdx, 


the diffusion energy 


-'diff 


{u, v):= (u log u + V log u) dx if m = 1 and 

jRd 


£dm{u,v) := 


1 


m — 1 

and the external potential energy 


(u™ + v^) dx if m > 1 


(2.4) 


£e^tiu,v):= / {uU + vV)dx. 


(2.5) 


Note that, with our assumptions, are finite for all {u,v) (z K x K. For 

"(/; = (—A)“^(ti — v) = G * {u — v) we define now the coupling energy 


£cp\{u,v) ■=^ [ iVV’pdx, 

^ jRd- 


( 2 . 6 ) 


which is the energy of the self-induced electric potential. Note that, at least formally, 
iVV'l^dx = / A'ij;)'ip dx = I {u — v) G * {u — v) dx 


[ |VV^|2dx= [ 
jRd- 


[u — v\{x)G{x — y)[u — v\{y) dxdy 


(2.7) 


falls into the category of interaction energies 



p{x)K{x,y)p{y)dxdy 


treated in |35| . To sum up, the total free energy £ = £dis + £ext + ^^cpl is given by 

(m 

For given probability measures p,iy ^ P(]R'^) we denote the squared (quadratic) 
Wasserstein distance by 


d^{p,u) 


inf 

7er(p,i/) 


X - yfdj{x,y), 


where r(/r, u) C x M^) is the set of admissible joint distributions with x and y 

marginals p,iy respectively. We recall from |35| that {'P,dw) is a metric space and 



that dw metrizes the weak convergence of measures. When fi € is moreove? 
absolutely continuous with respect to the Lebesgue measure d/x(x) <C dx, the square 
Wasserstein distance can also be computed by the Benamou-Brenier theorem |5] as 
already noted. 

Theorem 1 (Existence of optimal maps, |3S])- Let /i € and v £ ^^^(IR'^). 

There exists a unique optimal transport map T = Vy? £ L^(]R'^;d/r) for some convex 
function ip such that 

n = T^lx: V/£ Cc(M'^), / f{y)du{y)=f f oT{x)dfi{x) 

jRd JRd. 

and 

d‘^{p,,n)= [ \x -T{x)\‘^dy.{x). 


Our interest here is a system so we endow x with the natural product 

distance 

d‘^{z, z') = d^iu, u') + d^{v, v') 

for all z = (tt,u) and z' = {u\v') in V x V. It is well known |30[ [35] that {'P,dw) 
enjoys a natural differential structure defined by means of continuity equations, so 
that our product space also has the same differential structure. This permits us to 
differentiate real-valued functions J- on the product space, and defines the corre¬ 
sponding Wasserstein gradient by the chain rule ^iF{zt) = grad^y We 

show now that (HU is really the gradient flow 

z(«) = (“{(}). 


In terms of ordinary calculus of variations, we recall that this is achieved by variation 
of domain or, in fluid dynamical terms, by Lagrangian variations, |16j . To this end, 
let us split for convenience the coupling Poisson equation —Aijj = u — v as 


Ip = ipu — with 


f -Afju = u ^ il)u = G*u, 
y —Afty = V 'ipy = G * V. 


Formally differentiating the diffusive energy (|2.4I) with respect to u (resp. v), a by 
now classical computation [301 [35] leads to the Au^ (resp. Au”*) term in (|l.ll) . 
Similarly, differentiating the external energy (12.5p with respect to u and v classically 
gives rise to V • {uS/U) and V • {v'VV) in (11.11) . In order to differentiate the coupling 
term we first use the formal integration by parts (|2.7p and then exploit the symmetry 
G{x — y) = G{y — x) to expand 


^cpl(^) v') 


1 

2 

1 

2 



v]{x)G{x 


y)[u-v\{y) dx dy 


\u{G *u) + v{G * x)] dx 


u{x)G{x 


y)v{y)dx dy. 


Differentiating with respect to u, it is well known 135] that the first integral gives 
the corresponding V • (riV(G * u)) -|- 0 = V • [uS/ipu) term. Rewriting the remaining 
cross term 


IL 


uGv dx 




u{G * v) dx 



dx. 


and noting that pty is independent of u, it is again well known that this term gives 
rise to —V • (uS/'ipy). Summing up we obtain V • (uV {ipu — ipy)) = V • (uS/fj) as in 


fhe first equation of (11.11) . Similarly differentiating with respect to v we obtain the 
— V • {v’S/'i/j) term appearing in the second component. 

Though very general notions of solutions related to Energy Dissipation Equality 
(EDE) or Evolution Variational Inequality (EVI) can be used for abstract gradient 
flows in metric spaces, mm, we use the more direct framework, introducing some 
features later for the implementation of the flow-interchange method. 


Definition 2.1. A pair u,v : (0, oo) x is a global weak solution if u,v G 

C([0, oo)]V), u{t),v{t) tt°, in {V, dw) as t\0, Vu™, Vv^, uVU, vW, uVif: 
and vS/ijj G L^(0, T; for all T > 0, and for any fixed ip G 


_d 

dt 



dt 


v(t, x)ip dx 


[ (Vu™(t,x),V<y9)dx- [ u{t,x){VU,Vg?)dx 

jRd Jffid 

/ u{t,x){V'tjj{t,x),Vip) dx, 

[ {Vv^{t,x),Vp>)dx - [ v{t,x){VV,Vp})dx 
/ v{t,x){V'if{t,x),Vip) dx 


( 2 . 8 ) 


(2.9) 


hold in the sense of distributions 'D'{0,oo) with f;{t,x) = G*[u — v]{t,x) a.e. {t,x) G 
(0, oo) X M'^. 


We observe that oo); L^(M'^)) could be replaced by L)g^((0, oo) x M'^) in the 

above definition, which is enough for all the integrals in (I2.8I) - (I2.9D to make sense. In 
any case the weak solutions constructed here would still enjoy strong regularity in 
the end, and our choice of including the regularity in the definition of weak solutions 
is purely practical. 

In this setting our main result is 


Theorem 2 (Existence of solutions for m > I). Fix m > 1 and initial data u^,v^ 
as in ([ 23 ]). Then there exists a global weak solution (u, v) with 

u,v€L’^ (0, oo; n (I + |xp)dx)) 

S7i> G L°°(0,oo;L2(M'=')) 

12 mi 

771 m o / ^ 

u^,v^ e 

u,ve L°°(0,r;LP(M'^)), Vp G [l,2d/{d + 1)], 

for all T > 0, and 

£{u{t),v{t)) < £{u^,v^) for a.e. t > 0. (2-II) 

If we further assume that u^,v^ G L'^fR'^) for some p G [I, oo], then Vr > 0, 

SUP] (||^(0IIlp(R‘*) + 11^(0 IIlp(R‘^)) ^ ^||u°||i;,P(Kd) + ||'y°||LP(Rd)) , (2.12) 

with 

A = max III A[/||j;^oo(-Rd), ||AV||^oo(-Rd)}. (2.13) 


In the case of linear diffusion we have similarly 


Theorem 3 (Existence of solutions for m = 1). The conclusions of Theoreml^ holt 
form =1 if we replace u,v ^ L°°(0, oo; in (|2.10p by 

u,v € L°°{0, oo; logL^(M'^)). 


We would like to stress again that estimate (|2.12p holds for p = oo, and was not 
known in the whole space as far as we can tell. Since we interpret as a gradient 
flow one could expect energy monotonicity £{t) f. This would immediately follow 
from (j2.1ip and uniqueness of solutions. Unfortunately due to the lack of regularity 
and displacement convexity we were not able to prove uniqueness within the above 
class of weak solutions, and therefore we only retrieve an energy upper bound. 

It is worth mentioning that the gradient flow structure of the PNP system and 
the above theorems are also valid in the bounded domain case, with some mild 
assumptions on the boundary and minor modifications of the proofs. Suppose U C 
is a smooth bounded and convex domain, and consider the physically relevant 
boundary condition, that is, the no flux boundary condition 


du dv dip ^ 

dv dv dv 


on 


dTl, 


J 

Jdn 


ipdx = Q, 


where v is the unit outward normal on dTl. We also assume the external potentials 
satisfy 

dU dV 
— = — = 0 on 
dv dv 

By m we know that the electrostatic potential can be represented as 


ip{x) = / N{x,y){u-v){y)dy, 
Jn 


Vx G U. 


Here the (singular) kernel N{x,y) = N(y,x) serves as a counterpart of the Green’s 
function G{x—y) in for the Newton potential. Then we may argue in a similar way 
that the PNP system formally possesses the gradient flow structure, and the existence 
theorems can be proved in a similar but somewhat technically easier manner. 


3. Study of the energy functionals 

In this section we study various properties of the two relevant energy functionals, 
namely the total free energy £ and the functional (13.3p used in the JKO minimizing 
scheme. As already mentioned in the introduction, we use the flow interchange 
technique to establish some improved regularity for the minimizers, which will turn 
out to be crucial in the next sections. 

For further use we recall here a particular case of the celebrated Hardy-Littlewood- 
Sobolev (HLS) inequality: 

Lemma 3.1 (Hardy-Littlewood-Sobolev, |23l I33|l. In dimension d > 

LP(M‘^) and ^ = G * w. Then 

(1) //l<p< d/2 there is a G = G{p, d) such that 

11*^11 -JP-, — C'llw'llLP(Rd), 

while if p = 1 there is a G = G{d) such that 

||d>|| d < 

(2) If 1 < p < d there is a C = Cp^a such that 

||Vd>|| dp ^ Cllrull 

i3^(Rd) ^ > 


3 let w G 

(HLS-1) 

(HLS-2) 


(HLS-3) 








since |G(3:)| = and |VG(x)| = -pprr this is a particular case of well known 

fractional integration results for the Riesz potential /«/ = | |j_„ * / with a = 2,1, 

and we refer to |23l [33] for details. Here denotes the weak-L^ space and 

coincides with the usual Lorentz space 

As an immediate consequence we have the following integration by parts formula: 

Proposition 3.1. Let d> 3 and w £ Then $ = (-A)-iu; = G*w 

satisfies $ G G and 

/ |Vd>pdj; = / d>r(;dx= // w{x)G{x — y)w{y)dx dy. (3.1) 

We shall use this later on with w = u — v in order to control = G * {u — v). 

Proof. Taking p = -^ £ i'i-,d/2) in (IHTH-1 yHTH-311 we see that $ G L 2 d/(d- 2 ) 
and V$ G Since {2d/{d + 2))' = 2d/{d — 2) all the integrals in (|3.1I) are 

absolutely convergent and the last equality holds by Fubini’s theorem. In order to 
retrieve the first equality we use approximation: if G C)?°(]R'^) converges to w in 
j 2 d/{d+ 2 )then by the HLS lemma —>■ <I> in and ^ V$ in 

L‘^{Rfi). Since ()3.ip holds for smooth Wn £ with — = Wn we conclude by 

letting n ^ oo. □ 


Back to our energy functional, we begin with a fairly standard type of result 


Proposition 3.2 (Energy lower bound). Let m > 1 and fC as in (|2.2I1 . The total 
free energy £ is a proper functional on 1C x 1C and 

inf £(u, v) > —oo. 

Kxic 

Moreover we have in every sub-levelset {£{u,v) < R} that 

(i) gradient eontrol: ||VV'||£, 2 (]jjd) < G 
(a) no eoncentration: if m > 1 then 

[ {u^ + v^)dx <G, 

dRd 

while if m = 1 then 


(u| log u\ + uj log u|) dx < G. 


(in) mass confinement: |xp(ri + v) dx < G, 

for some G > 0 depending on R> d, the confining potentials, and m. 

Proof. Choosing u,v smooth and compactly supported it is clear that £{u,v) < oo 
so £ is proper. 

m > 1 : (i)-(ii) immediately hold because each term in (II.3h is nonnegative, (iii) 
then follows by fj^a(uU + vV) dx < £(u,v) together with (12.ip . 

m = 1 : if rri2(/9) = |x|^/9dx denotes the second moment let us first recall [16] the 

Carleman estimate 

R(p) >- [ P(log P)~ dx > -G{1 + m2{p)T,p £ V, 

jRd 


(3.2) 









for some C > 0 and a G (0,1) depending on the dimension d only. By (|2.ip we ha"^^ 
Sext{u,v) > Ci{m 2 {u) + m 2 (r')), whence 

C [m2(ri) + m2(u) - (1 + m2(rt))" - (1 + m2(?;))"] 

< 1-L{u) + T-L{v) + [ {uU + vV) dx < £{u, v) < C. 

Hence the second moments are bounded as in (hi). Then (i) and (ii) come immedi¬ 
ately from (hi) and (|3.2p . □ 

For fixed 2 ;* = (tt*,n*) ^ V x V, and given time step h > 0 we set 

I'hiz) :=-^d‘^{z,z^) + £{z), z = {u,v) e 1C X 1C. (3.3) 

In order to dehne later a discrete sequence of approximate solutions using the JKO 
minimizing scheme, we collect here some properties of and preliminary results. 

Proposition 3.3 (Existence of minimizers). Fix h> 0, and z^, = (tt*,u*) G P x P. 
Then IFh admits a unique minimizer z = {u,v) (z K. x 1C. 

Proof. By Proposition 13.21 Th is bounded from below on 1C x 1C, hence there is 
a minimizing sequence Zk = {uk,Vk) satisfying (i)-(iii) and {uk,Vk}k are tight and 
uniformly integrable. By the Dunford-Pettis Theorem one may extract a subsequence 
such that 

Uk ^ u and Vk ^ v in L^(]R‘^), 

and standard truncation arguments together with the uniform bounds on the sec¬ 
ond moments ensure that u,v £ V. The weak lower semi-continuity (l.s.c.) of 
the squared Wasserstein distance, diffusive and potential energies are standard, in 
particular u,v £ 1C. We prove in the appendix, Proposition 16.11 that the Dirichlet 
energy is lower semicontinuous with respect to weak L^(]R'^) convergence. Because 
Uk — Vk ^ u — V in L^{W^) we conclude here that £cn\{u,v) < liminf £lcpl(^fc) 

thus u,v is a minimizer. Finally, the uniqueness result comes from the fact that the 
admissible set /C x /C is convex w.r.t. linear interpolation zq = {1 — 6)zq + Qz\ and 
that the total free energy is jointly strictly convex in {u,v). □ 

We remark that the squared distance term left aside in (13.3p . the same line of 
argument would readily give existence of a global minimizer of the total free energy 
£, which would result in the end in a least energy stationary weak solution u,v to 
m- Since our gradient flow system is driven by £ one could expect long time 
convergence u{t),v{t) u,v when t —)• 00 together with some convergence rates. 
However, the lack of displacement convexity prevents here from applying standard 
techniques mull [30] and this is beyond the scope of this paper. We refer to El El 
for related results on similar PNP models. 

In the next section we shall derive the discrete Euler-Lagrange equations satisfied 
by the minimizers, which requires integration by parts as in ()3.ip . However, at this 
stage the minimizers only lie in L™'(M'^) if m > 1 and log L^(M'^) if m = 1, and this 
manipulation is not justified. The discrete Euler-Lagrange equations are necessary 
to pass to the limit as the time step /i —> 0 and to thereby obtain a solution to the 
PNP system. The remainder of this section is devoted to improving the regularity 
of the minimizers of the discrete functional. 

The argument is based on the flow interchange technique of Matthes, McCann, and 
Savare, [25], as implemented by Blanchet and Laurengot m, as well as Laurengot 
and Matioc m and El . The idea of the flow interchange technique is that a known 


^'adient flow is sufficiently close, generally first order close, to the one under study 
so that it may be used as an approximation with controllable error. We need to 
use this method twice, first to propagate the regularity of the minimizers and then 
to establish some smoothness of their spatial gradients. The characterization of 
gradient flow that is useful here is the so called Evolution Variational Inequality 
(EVI) for a functional T. Using the notation to follow, a flow {u{t)) is a gradient 
flow in the EVI sense, mm, provided that 

re) + V'(n(t)) < J-'(rc) for all w G and a.e. t > 0. (3.4) 

Displacement convexity and the other detailed requirements for (|3.4I) to hold are 
discussed in the references just cited. Eor our purposes we note that (13.41) is valid 
for 

(1) solutions of dtu = Au, the heat equation, with V = "H, the Boltzmann 
entropy 

■H(fl) = / ulogudx and (3-5) 

(2) solutions of dtu = A (u^), the porous medium flow, with = £p, 1 < p < oo, 
given by the functional 

Spin) := -- [ iFdx. (3-6) 

P — ^ jR<i 


Proposition 3.4 (Discrete propagation of estimates). Let m > 1, X as in (|2.13p . 
and further assume that tt*, u* G }CriL^{W^) for some p G (1, oo). If0<h< ho{p) = 
X(^p_i'j then the minimizer {u, v) from Propositionsatisfies 


wr 


LP 


+ T 


LP 


< 


1 


1 — X{p — l)h 


(ii-.il 


+ P'* 


IP 

Ilp 


(3.7) 


In this first use of the flow interchange, we simply use the solution of (|3.8I) below 
as variations in the minimum principle. Note that at this point the time step h must 
be taken small in terms of p for the minimizing problem to “see” the estimate. As a 
consequence there is no hope to retrieve an estimate at the discrete level for 

fixed h directly from the limit p ^ oo in (|3.7p . since h < hQ{p) would require /i —0. 
However, u, v will be retrieved as some limit when /i —>■ 0, so one can actually take p 
arbitrarily large and the weak solutions will ultimately satisfy such an L°° estimate. 
See the proof of Theorem [2] at the end of Section O for details. 

Proof. Eor fixed p G (1, oo) and u*, u* G JCflL^ consider the auxiliary Porous Media 
flows u{t),v{t) defined by 

dtu = A (fl^) in (0, oo) x u\t=o = u in 

dtv = A (v^) in (0, oo) x v\t=o = u in ^ ^ 

By standard results for the PME |34) we know that (i) these Cauchy problems 

are well posed and u,v € C([0, oo); (M*^)) remain probability measures, (ii) by 

pi _ j^oo snioothing u(t),v(t) G for all t > 0, and (hi) the second moments 

remain finite. As a consequence u{t),v{t) G JC are admissible for any t > 0, and by 

Proposition 13.11 it will be no issue to integrate by parts in the coupling term. 

Step 1: dissipation of the internal energy. We first claim that 

d /- \ 

—£dm{u,v) < 0 , 


t > 0 : 


(3.9) 
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and we distinguish cases, m = 1 being the most involved, 
m > 1 : By usual properties |34) of the PME all the norms are non-increasing 

along the PME flow, in particular £dm{u,v) = (ll«lir-‘(Md) + ll^llL-^(Md)) is non¬ 
increasing in time. 

m = 1 : Assume first that u is smooth and positive. Then by standard properties of 
the PME flow so is u{t) for later times and, therefore, we get 

f ulogudx) = [ (logu)Au^dx = — - [ \Vu^^‘^\‘^dx < 0 (3.10) 

di jRd. / Jjjd p J^d 

for t > 0. Thus 1 1 —)■ T-L{u{t)) is nonincreasing. 

If u is not smooth and positive we first regularize it by running the heat equation 
for small times 

Uk=Ti/k*u u, 

k^oo 

where P^ is the usual heat kernel at time s. Since the heat equation is the "H-gradient 
flow we have in particular T-L{uk) < 'H{u), and of course Uk is positive and smooth. 
Denoting by Uk{t) the solution of the corresponding PME-flow dtUk = starting 
from Uk, then the previous computation (j3.10p shows that 

vt>o: niukit)) ^niukio)) ^niu). 

Since Ufc —)• u in we get by standard contractivity of the PME that 

Uk{t) —>■ u{t) in uniformly in t > 0 when k —>■ oo, in particular weakly in 

By lower semi-continuity of Ti with respect to weak convergence we 
conclude that 


Vt > 0 : 'H{u{t)) < liminf'H(ufc(t)) < 'H(u). 

k^oo 

Finally, by uniqueness of solutions to the PME flow with u{0) = u we conclude that 
t I—> is monotone nonincreasing, and similarly arguing for v entails (13.9h as 

claimed. 

Step 2: the remaining terms. Arguing by approximation |34| the potential 
energy is easily controlled for t > 0 as 

^£e^t{u,v) = ^ f {uU + vV) dx = f [{AuP)U + {AvP)V] dx 
at at jRd J^d 

= [ \tFAU + v^AV]dx<\ [ {iF + v^)dx. (3.11) 

jRd JjRd 

For the coupling term, let ^(t) = G * {u — v){t) and observe that for t > 0 we have 
dt'ip = clt[(—A)“^(u — f})] = (—A)“^[(9t(u — u)] = —{uP — yP). Since € 

for t > 0 we can legitimately integrate by parts 

4£,p,(i.ii) = 4(^JV^pdx 

= f {—A'tp)dt'il’dx = — f {u — v).{u^ — v^) dx < 0. (3.12) 

JRd 

Note that, due to ||V'0||L2(]Kd) = ||V(—A)~^(u — u)|| 2 , 2 (]gd) ^ \\u — 'i}||j:^-i(rod',, this is 
the well-known contraction property of the PME flow, see 





As for the Wasserstein term, note that u, v are respective gradient flows of the 
functional £p, (|3.6I) . so from (|3.4I) . 


+ dw{v, 

< 7 —[ « - uP) dx + - —3— [ (yP - yP) dx 
[P - l)/l jRd [P - 1)^ jRd 


(3.13) 


Step 3: dissipation inequality. Gathering (13.91) . (|3.1ip . (13.121) . and (13.131) . we 
get the total dissipation inequality 


V(t) := 

+ a/ 


/ (u^ + v^)dx - - —1— [ (u^ + v^) dx 
jRrf [P - Ij/i Jmd. 

(yP + yP) dx =: A(t) 


for small t > 0. Because (t((0), i;(0)) = (u,v) is a minimizer we must have D(t) > 0 
at least for a time sequence \ 0, otherwise {u{t),v{t)) would be a strictly better 
competitor for small t > 0. If 0 < h < ho = x(p-i) have 1 — X{p — l)h > 0 and 
A{tn) > TX{tn) > 0 can be rearranged as 

[ [u^{tn) + vP{tn)] dx < ---— [ {yP + yP) dx. (3.14) 

jRd 1 - A(p - l)h Js^d 

Our statement follows by finally letting 0 in (|3.14p . recalling that (u(t), i}(t)) —>■ 
{y,y) in L^{W^) when t —)• 0. □ 


We shall also need a further regularity result for the gradient of (u, y). The use of 
the flow interchange in this estimate is very similar to its use in |10j for the critical 
parabolic-parabolic Keller-Segel model. 


Proposition 3.5 (Discrete gradient estimate). For m > 1, d > 3, and any h > 0, 
fix z* = (u*,x*) € V X V and let z = {y,y) IC x 1C be the unique minimizer from 
Proposition \S.tA Then 


< c 


h 

for some G > 0 independent of h > 0 and z^,. 


^ ^ n{u,)-n{u) ^ n{y,)-n{y) 


(3.15) 


Proof. We use a second flow interchange with u{t),v{t) now defined by 

dtu — Ay = 0 in (0, oo) x y\t=o = u in (3.16) 

and 

dtv — Av = 0 in (0, oo) x R'^, v\t=Q = x in R'^. (3.17) 

Step 1: dissipation ineqnality. We first note that classical properties of the 
heat equation and n(0),x(0) G 1C guarantee u{t),v{t) G ^ for all t > 0. Let := 
G * [y — v). Then it is easy to check that dtfi> = Afi as well. Since the pair (rt, x) 
is a minimizer and has finite energy we have in particular V'i/’(0) = V'i/’ G L^(R'^), 
whence by standard properties of the heat equation V'i/((t) G L^(R'^) and 

<0 

for all t > 0. Since y, v are positive and smooth for all t > 0 we may differentiate 
and integrate by parts as 
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m > 1 : 




m 


m — 


- [ {u^-^Au + v"^-^Av)dx+ [ {UAu + VAv)dx 

1 jRd J]^d 


< —— [ + |Vi}'2'p')dx + [ (uAU + vAV)dx 

'm Jud\ ) J^d 

~ ~m j d ^ P + 1^^ ^ P) dx + ||AC/||^oo(]Kd) + ||AV||^oo(Rd) 


(3.18) 


and similarly 

m = 1 : 

—£{u,v) < ^iVtia |2 -I- |V'i} 2 dx + ||A[/||^oo(Rd) + ||AF||j;^oo(Kd). (3.19) 

Here (12.11) is used in the third inequality of both (I3.18P and (|3.19p . 

For m > 1, the above term ^/jRd | Vp” pdx in p3.18p (dissipation of £m{p) = 

Fiti /iRd P"* dx along dtp = Ap) corresponds to the usual Fisher information —4 |Vp 2 \‘^dx 
^'H(p) in the linear diffusion case p3.19p . and enjoys a formal continuity when 
m 1. Comparing with (I3.10p for p = m we also see that the dissipation of Ti 
along the £m-^ow equals the dissipation of £m along the ^-flow, which is in fact the 
cornerstone of this flow interchange technique. 

Let 

m := 4 (||Vi¥ + II . 

Integrating (j3.18p or (j3.19p from 0 to t > 0 we get in both cases 

£{u{t),v{t)) — £{u^,v^) < 2Xt — f 'D{s)ds, (3.20) 

Jo 

with A defined in (j2.13p . Because ()3.16p - (j3.17p are respective "H-gradient flows, we 
again appeal to (13.4p to obtain 


< L^(tt*) — + T-L{v^) — 


Integrating again and using the monotonicity of s \ 'H{u{s)) and s T-L{v{s)) along 
the flow with 5(0) = z gives 

^ [d‘^{z{t), z^) - d^{z, z*)] < [H{u^) - 'H{u{t)) + 'H(x*) - . 

Since 2 ; is a minimizer we have by (13.3h and (I3.20p that for small t > 0 

0 < d^h{z{t)) - Thiz) 

<— V.{u{t))+T-L{v^:)+2Xt — / D{s)ds, 
d Jo 

which we reformulate as 

1 f‘ V{s) ds< 2 X + «(“■) - «("(«» + «(”•)-«("(*)). ( 3 . 21 ) 

t Jo h h 




























^ep 2: the limit t —> 0. If 

/ , 1 ^ ^ 1 ^ 

Diit^x) ■.= - / ft 2 (s,xjds, D 2 [t,x) := - / i; 2 (s,xjds, 

^ Jo ^ Jo 

we first note that u,v ^ C([0, 00 ); L™'(M'^)) as solutions of the heat equation with 
initial data in so that Di,D 2 G C([0, 00 ); L^(M'^)). As a consequence 

Di{t) —)■ -Di(O) = u~ and D 2 {t) —>• 1 ) 2 ( 0 ) = in when t I 0. By (|3.2ip 

we find that VDi(t) and SID 2 {t) are bounded in and converge at least in 

to V(u^) and V {y ^) when t —>■ 0. Consequently, V {u '^), V) G 
and our statement follows from 

II ^ )|lL2(Rd) + ||V('f^ ^ )|L2(Rd) 

< li^nf ||VDi(t)||5^2(^,) + li^nf ||VZ)2(t)||i2(R.) 

Tfl 1 

< — liminf- / V(s)ds 

4 t^O t Jq 

^ ^ ^ nM-nju) ^ n{v.)-n{v) \ 

“ \ h h J ' 

The last inequality comes from the limit t \ 0 in (I3.21h with strong convergence 
{u{t),v{t)) —>■ {u,v), e.g. in n and suitable continuity of p 1 -^ TJ{p)- D 

4. Minimizing scheme and discrete estimates 

In this section, we shall construct a family of time-discrete approximate solutions 
using the JKO method, also known as the variational minimizing movement scheme. 
A priori estimates for the set of discrete solutions are necessary to allow us to deduce 
the existence of a time-continuous limit curve. 

Fix an initial datum = (u^,v^) as in (12.3h and some time step h > 0. Setting 
Proposition 13.31 allows us to dehne a sequence ^ K, x K, 

recursively as 

:= the unique minimizer 2 : of Fh with 2 * = 2 ^”^ = 
and a corresponding piecewise-constant interpolation t G [0, 00 ) i-A Zh{t) as 

Zh{t) = z^ff^ for nh <t < {n + l)h. 

The rest of this section is devoted to collecting the suitable a priori estimates on Zh 
suitable to pass to the limit in /i \ 0. 

It is now standard to get the discrete energy monotonicity as in |16) that 

V„>0: S(4”+‘>) < «(4”’) 

inasmuch as z^'^ is a competitor in the search for At the continuous level 

this reads as 

Sizh{t 2 )) < £{zh{ti)) < £{z^) for all 0 < ti < t 2 - 
Proposition 4.1. The total square distance and approximate Holder estimates 

(4.1) 

n>0 

VO < ti < t 2 , d{zh{ti),Zh{t 2 )) < C\t 2 -ti - 4 /l |2 , 
hold for some C > 0 independent of h > 0. 


(4.2) 








Proof. Note that since u^,v^ G we have in particular — 

j^ 2 d/(d+ 2 ) gy Proposition 13.1 1 we have Vijj^ = VG*{u^ —v^) G and u^,v^ 

has therefore finite energy. We also recall from Proposition 13.21 that > — oo, 

hence the right-hand side in (14.ip is finite (and of course independent of /i > 0). The 
rest of the argument is by now very classical and we refer to |16) . □ 

Proposition 4.2. The piecewise constant interpolation satisfies 


m > 1 : sup f (u^{t) + v^{t)) dx < C 

t>0 jRd. 

m = l:sup / (n/j(t)|logUfe(t)|-bu/i(t)|logUfe(t)|) dx < (7 
t>0 jRd 


(4.3) 


and 


sup [ 1x1“^ {uh{t) + Vh{t)) dx < C (4.4) 

t>0 JRd 


uniformly in h > 0. 

Proof. By energy monotonicity we have sup£{u^f^\v^'^^ < £{u^,v^) < oo, which by 


n>0 


Proposition 13. 2l bounds the internal energy and the second moments uniformly in h, n 
for the discrete sequence. This property extends to the interpolation Uh{t), Vh{t). □ 


In addition to the uniform control in Proposition 14.21 we also have 

Proposition 4.3 (Continuous estimate). In addition to (12.31) assume that the 
initial data u^,v^ G LP(M'^) for some p G (l,oo), and let X as in (|2.13p . Then for 
h < ho(p) = sufficiently small we have 

Vt > 0 : |h/i(t)|Lp(Kd) + lk/»(0|Lp(Kd) < Ce^^(\\u°\\LP(ud) + ll^^°llLP(Erf)) 
for some C > 0 independent of t, p, h, and the initial data. 

Proof. Fix any t > 0, let /c = [t/h\, and recall that Uh{t) = u'^'^. By induction we 
immediately get from Proposition 13.41 

For small h > 0 this easily gives 

II'“/i(^)IIlp(R‘^) + ll''^/i(^)llLP(Rrf) ^ P *(||'«°||lp(R'^) + II^°IIlp(R'^)) 

for some universal C > 0. Since e p < the proof is complete. □ 

Proposition 4.4 (Approximate Euler-Lagrange equations). Fix m > 1. Let 
and Vr*'"') be the optimal transport maps 






in Brenier’s Theorem m and = G * (uj^^ — Then for any vector-field 


C G we have that 


-Id,Co Vg(”Vjr^dx 


[ (^^Jr^^^)”'divCdx- [ u^f^^^''{VUX)dx - j Q) dx, 

jRd JRrf jRd 


(4.5) 


and 


1 f (Vr(") - Id, C o Vr^^^)vl^'>dx 


h 


[ div Cdx — f C) dx + f dx. 

jRd. jRd jRd 


(4.6) 


Proof. In order to simplify notations, we write below rt* = u^i^\ u = x* = 

V = and fi = arbitrary vector-field 

C G C^(M'^,M'^). For e G [—(1,(1], let 4>£(x) be the associated e-flow (i-e d^e/de = 
^(^e) and <ho = Id), and let us consider the perturbation (of domain) 

Ue-={^e)#U, Ze-.= {Ue,v). 

Since z\e=Q = 2 is a minimizer, computing the hrst variation ^ {^hi^e)) i;=q = 0 
will classically give (|4.5I) . Similarly considering = (‘hg)^^ and = {u,vfi) will 
produce (14.6p . 

More precisely, differentiating the Wasserstein distance squared, the confining 
potential, and the diffusive energy are by now classical computations [2]. However, 
differentiating the coupling energy is quite delicate here: because we have to consider 
separate horizontal and vertical perturbations the nonscalar nature of the problem 
induces a loss of symmetry. Formally the result should follow from 

/ \Vfie?dx= / fie{u£- v)dx = // [ue - v\{x)G{x - y)[ue - v]{y) dxdy 


and the classical computations for interaction energies, see |35j . But because we 
consider two components independently it might happen that ^ L^(M'^) even 
though VV’ G L^(M'^), and the above integration by parts might not be legitimate. 
Moreover since VG is more singular than G itself, differentiating with respect to 
e requires some extra regularity. This can actually be made rigorous using the 
propagation of the initial regularity as follows. Since the initial datum vP G 
j^ 2 d/(d-ei) and the time step is small enough, we have by Proposition 13.41 that 

u,v G L^{W^) D and in particular u,U£,v G (M'^). Using 

Proposition 13.II we can therefore integrate by parts and expand with = (4*£)^u 


\Vfie\‘^dx 


u{x)G{^fix) - ^e{y))u{y) dxdy 


- 2 


// 




u(x)G(‘he(x) — y)v{y) dxdy -\- terms independent of e. 


where the last equality follows by dehnition of the pushforward In 

order to differentiate under the integral sign we only need L^(M'^ x R*^) bounds such 


that 


19 


u(x) (VG(^s(x) - ^e(y)), C o -Co •^e(y)} u{y) dxdy < C, 

__ xR'* 

[[ u{x) {VG{^e{x) -y)x° ^e{x)) v{y) dxdy < G, 

uniformly as e ^ 0. Because <!>£ is close to Id for small e, C ^ and 

this simply amounts to controlling 

11 ^ c*’ 

JjRdx-Rd \x-yr ^ 

11 u{x)- - K^v{y) dxdy < G, 

JjRdy,^d \x-y\'^ ^ 

which is valid by (IHLS-3P with p = 2d/{d + 1) and u, u G As a 

consequence we can legitimately compute with <ho = Id 
d 


de 


([ IViA.pdx) =[[ 

^ jR'i ^ e=0 J JR' 


u{x){VG{x - y),C{x) - C{y))uiy) dxdy 

xR'* 

- 2 //' u{x){VG{x - y),C{x))v{y)dxdy. 

J JR'^xW^ 

Exploiting the symmetry VG(x — y) = —'S/G{y — x) we hnally get 
l^dx 


U\i 

de 2 jjjd 


u{x){VG{x — y), C{x))[u — v]{y) dxdy 
= / u(V'0,C)dx 


as in our claim, and the proof is complete. □ 

The above restriction at initial data u^,v^ G (1^^)) which then is inherited 

by the solutions to later times, is technically essential in order to differentiate under 
the integral sign with respect to e-perturbations and retrieve the discrete Euler- 
Lagrange equations. Actually this restriction is not purely technical: in (HH) it 
seems natural to require the terms uVip, vVijj to be at least in at time t = 0. 

If u^,v^ are both in for some p then the integrability for coming from 

(IHLS-3P is Vi/j G , which is optimal since HLS inequalities are. Solving 

for p' = gives exactly the sharp p = 2d/{d + 1) exponent. Technically speaking 
we had to assume initial regularity with slightly better but arbitrarily close 

ro > 2d/{d + !)■ This is needed for technical compactness issues, arising later on 
when we take the limit h ^ 0 to retrieve the weak solution (tt, v) = lim(rt/i, Vh)- 

In addition to being an approximate solution in the sense of the previous Propo¬ 
sition, the interpolation {uh,Vh) satisfies 


Corollary 4.1 (continuous gradient estimate). Fix m > 1. Then for allO < h <T, 
||V(n/i) ! + ||V(u/i) ! < C'(^ + 1) 2 , (4.7) 

and 

||V(u/i) ||i2(;j 2’;L1(R‘*)) IL2 (/i_T;L1(R‘*)) — 

for some constant G = G{vP,v^) > 0 independent of h. 













l^roof. For once the argument requires no distinction between m > 1 or m = 1. We 
only estimate the u component because the computations are identical for v. 

Since G we have 

m m 

V(u,r = 2<V«) gL1(M''). 

Recalling that is actually bounded in uniformly in f > 0 and h, 

clearly (14.Sp will follow from (|4.7I) and we only establish the latter. 

For fixed 0 < h < T let N = [T/h\, and recall that the interpolation Zh{t) is 
piecewise constant. Multiplying ()3.15p by /i > 0 and summing from n = 0 to n = N 
we obtain 

. 

Ih 

Af-1 


|V('U/i) ■ ||^2(Rd) 

riN+l)h 


< 


< 




L2(Rd) 


n=0 


L2(IRd) 


Af-1 


C (a+«( 4”>) - «(4”+‘>)+«(4”>) - w(4”+‘>)) 

n=0 

<c(t + n{u^) + n{v^) - - ^(4^^)). 


(4.9) 


By Proposition 14.21 the second moments m 2 




,m2 


) are bounded uniformly 


in f, h, n, hence by the Carleman estimate (|3.2I) we see that —Ti (D < c 

in ()4.9I) and the proof is complete. □ 


We observe that another possible way to retrieve better gradient regularity is to 
estimate 


/ div(CXdx 

jRd. 


+ 


div(C)urdx 


^ C'IIClll,P(Rd) 


for arbitrary vector-fields G M'^) in the Euler-Lagrange equations (14.5p - 

(14.6p . which would estimate by duality G see e.g. [29]. This 

approach would only improve the previous total variation estimate if p < 00 , so 
that {LP'= LP{'U.'^) and is dense in L^. Unfortunately we are here in 

a limiting situation where essentially UhV'iph,Vh'V'ipii G only, so this is not 

feasible. More precisely, our assumption u^,v^ G with tq > 2d/{d+l) in fact 

does give slightly better Uh'^iph G integrability through HLS inequalities 

(for some (5 > 0 depending on tq). But since tq could be arbitrarily close to the 
critical 2dj{d + 1) exponent, 5 > 0 is arbitrarily small and we shall refrain from 
taking this technical path. 


5. Convergence to a weak solution 

This section is devoted to the convergence of the previously approximated interpo¬ 
lating solution towards the final weak solution, {u,v) = llm{uh,Vh) in some suitable 

h —^0 

topology. Because of the quadratic interaction term and nonlinear diffusion if m > 1 
we will need the following strong convergence 

Theorem 4 (Strong convergence). There is a discrete subsequence, still denoted 
h\,0, and functions u, v such that 

Uh{t,x) ^ u{t,x) and Vh{t,x) ^ v{t,x) a.e. in (0, 00 ) x (5.1) 















and 
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VI < p < oo, 1 < O'< ro : Uh,Vh ^ u,v m oo); (5.2) 

where tq > max{m, 2(i/(d + 1)} is the initial integrability as in (|2.3p . 

Observe in particular that q = m and q = 2d/{d + 1) are allowed in (|5.2p . which 
will be crucial in order to pass to the limit in the Euler-Lagrange equations later on. 
Roughly speaking, L^{W^) integrability suffices to guarantee convergence for 

all q < r. Unfortunately q = max{m, 2d/{d + 1)} is a borderline case that we could 
not treat, and this is why we needed to assume initial integrability for some 

slightly better but arbitrarily close tq > max{m, 2d/{d + 1)}. 

Strategy of proof : We will use a compactness criterion in Bochner spaces from |32) 
that involves 

(i) boundedness in L^(0,T;X) for some strong X topology, 

(ii) compactness in LP{0,T]Y) for a weaker Y space, 

(hi) a target intermediate U’{0,T; B) space with embeddings X CC B CY. 

Note that X CC B will be achieved by space difference quotients and that the 
key ingredient to obtain boundedness of {uh,Vh}h in the strong X topology is the 
gradient estimate from Corollary 14.11 Compactness in LP{0,T;Y) will be ensured 
by an approximated time equi-continuity in some suitable (R*^) space. 


We first collect some technical results and then establish Theorem [H To begin 
with, for m > 1, we set 


1 777 — 1 

g^:=l + —= 1 + — G(l,m). 


m' 


m 


Then there exists 0^ G (0,1) satisfying 


and we let 


- — (1 “ dm)— + dm —; 

Qm I m 


2m 

Pm-= — > 1- 


dn 


If Te denotes the usual shift operator in space 


e G R“ : Tew{x) := w{x — e), 


./ 


(5.3) 

(5.4) 

(5.5) 


\x\^dx < oo 


we also define the weighted Nikolsk’ii spaces 

C 1 ^Tn 

Xm := iw € L'^"*(R ) : sup \\TeW — wWiqm \e\ ^ < oo 

eGK'* 

endowed with their natural Banach norms with dm!'^ < 1- By the Riesz-Frechet- 
Kolmogorov Theorem we have 

Xrr, CC L'''"(R'^). 


We note that the above choice for p = pm, Q = Qm^d = 6m is purely technical so we 
shall go as little as possible into details regarding their explicit values. 

For the case m = 1 one should similarly use 

Xi := \w G L^(R'^) : Xw G L^(R'^), [ \x\‘^\w\dx < oo] CC 

jRd. ) 

Since the related argument is fairly easy compared to the nonlinear case, in what 
follows we will omit the related proof and focus on the nonlinear diffusion case and 
henceforth we assume m > 1. Compactness in space will be ensured by 



Proposition 5.1 (Compactness in space). Let p,q,9,X as in (I5.3p -(f5]), and fix 
T > 0. Then for h > 0 small enough we have 

\\Uh\\LP{h,T-,Xm) + \\'^h\\LP{h,T-,Xm) ^ C't 

uniformly in h. 


Proof. For simplicity we write here p = pm, q = Qm, 9 = 9m, X = Xm- We first claim 
that 

\\uh{t)\\x < C{l + for t > h. (5.6) 

Indeed since q G (l,m), it follows immediately from (|4.3I) that 


ll^/i|lL°°(0,oo;L9(Kd)) < C. 

Using (|4.3I) . (14.4p . and Holder inequality we estimate with q = qm = ^ + ^/na' 

/ \x\'^u\{t) dx = / {uh{t)\x\^)^ Uh{t)dx 

< \\Uh{t)\\Lrn(^d) < C. 

Fixing e G and using the convexity inequality |a — < | |a|™ — |6|™'| for a, 6 > 0, 

we get 

\\TeUh{t) - \reUh{t) -M/i(t)rdx^" 

< ( [ \TeUhit) - U'^{t)\ dx) 

< \e\^/-\\VuT{t)\\%,y 

By (|5.4I) and \\TeUh{t) — n/i(t)||^i(Rd) < 2 we get by interpolation 

IliitiftW - “ii(*)||i,(ii,j) £ \\TM(t) - “(.(*)||l7(V)' - “'■WIIUpj) 

<2|e|^||V«r(()||f.(..,. 

thus (|5.6p holds as claimed. 

Taking now the IP{h,T) norm with p = pm = 2m/9 in (|5.6p and using Corol¬ 
lary [TT] finally leads to 

\\uh\\LP{h,T-,X) < C (t + \\Xui^ IL2(/i < C't- 

The estimate for the v component is again identical. □ 


Next, we turn to compactness in time in a weaker topology. We first have 

Proposition 5.2 (Time-equicontinuity in kF“®’'’^(]R‘^)). Let s,r > 0 be large enough 
so that 

Then 


Uhit2) - Uhitl)\\^_,y + \\vh{t2) - Vh{tl)\\^y_sy < C^/\t2 - tfi + h, 0 < H < t2, 

(5.7) 


for some C = Cg^r > 0 independent ofti,t 2 o.nd h. 



Proof. The argument is very similar to cni Lemma 13], beginning with the ca¥ 
culation of the approximate 1/2-Holder continuity of the sequence {uh{t)}. For 
indices 0 < n < n' let Vg denote the optimal map from to ^ so that 
Then 

[ {f(yq{x)) - ^{x))u]^\x)dx, f, 

jRd ^ ^ JRd 

Expanding the integrand on the right, 


Hence 


/ 


iix) - e(Vg(x)) = Ve(Vg(x)) ■ (x - Vg(x)) 

+ O(|x-Vg(x)|2||v2c||io 


4” ^ ? dx = / [? o Vg - ^] dx 

' JRd 

= [ (Vg - Id, VC o Vg)4”^ dx 

jRd 

+ 0(||V2c||ioo(R.)) / |Id-Vg|2uif)dx 

= [ (Vg - Id, VC o dx 

jR<i 


+ 0(||V2C||ioo(K.)dH^(ufUl"V) 

We further compute by Cauchy-Schwarz and Holder inequalities 


< 


(Vg-Id,VCo Vg)ui^”^ dx 
[ iVg-Idpuj/'^dx 


|VCo Vgpuj;'^ dx 


<dw{u’f^\u’'^''^)( I |VClX"Mx 




h ’ “h Hn/i IlL” 


< dw{ul \ ul ' ]\\u 


livci 


1 

2 || 2 




Given 0 < H < ^2 and Ni = , N2 = [t2/h \, from fid.lh and the Cauchy Schwarz 

inequality we get 


N2-1 

dw{Uh'^,Uh^) < ^ dw{u^h\u^h^^^) < C^/\t2 - ti\ -h h, 
n=Ni 


/ {uhit 2 ) -Uh{ti))fdx 
jRd- 


< c 


L2m' md)\/\t2 — Hi + ^ + I|V^C|Il° 



and then 








With our choice C l^h 2 m' p ^ 2,00 


and because h is small we finally obtain 


/ {uh{t2) - Uh{ti))^dx 

Jr<^ 


< C — ti\ + h + 

< C^/\t2 — til + 


Our statement follows by density of in 14^^’’'(M'^) and duality = 

□ 


We are now in position to prove the desired convergence when h ^ 0: 

Proof of Theorem\^ Once again we only establish the result for the u component. 
Fix any 0 < 6 < T and let q = q-m, 9 = 6m,P = Pm, X = Xm as in (I5.3l) -(f^. Taking 
s, r large enough such that 

CC 

is compact. By truncation, a standard duality argument then ensures that 

n L^(M'^)((1 + |x|2)dx) CC 

is also compact. By ProDosition l4.2l we see that there is a fixed W~^’'’^(M'^)-relatively 
compact set K such that Uh{t) G K for all t > 0 and small /i > 0. Therefore, we infer 
from Proposition 15.21 that a refined version of Arzela-Ascoli Theorem |2l Proposition 
3.3.1] can be applied to conclude that there exists u G C([0, T]; such 

that 

VtG[0,r], Uh{t)^u{t) in 

for some (discrete) subsequence /i \ 0, not relabeled here for simplicity. This point- 
wise convergence together with the uniform PI (1 -|- |xp)dx) bounds 

and Lebesgue’s Dominated Convergence Theorem therefore guarantee strong con¬ 
vergence 

u/, ^ u in TP(0,r; (5.8) 

By diagonal extraction we can moreover assume that u G C([0, 00 ); and that 

(j5.8p holds for all T > 0. 

Choosing s, r large enough we can further assume that 

A CC L'?(M'^) C W"’’''(M'^). 

We recall from Proposition 15.11 that {uh\h is bounded in LP{5,T] X), and by (15. 8p 
it is also relatively compact in L'^{6,T; (M'^)). By [321 Lemma 9] we conclude 

that {uh}h is relatively compact in the intermediate target space, i-e rt/i —>■ rt in 
L^(5, T; L'^(M'^)) for some subsequence. By a diagonal extraction we may assume 
that u is independent of 6, T, and Uh ^ u in 00 ; L'?(M'^)). Up to extraction of 

a further subsequence this classically implies the desired pointwise convergence a.e. 
(t, x) G (0, 00 ) X 

Let us turn now to the L](j^([0, 00 ); L'^(M'^)) convergence, and fix 1 < p < 00 
and 1 < (/ < ro as in our statement (we recall here that tq > max{m,2d/{d + 1)} 
is the initial integrability u^,v^ G L^°(M'^)). Once again we only focus on the u 
component. Thanks to the previous {t,x) a.e. convergence we shall apply Vitali’s 
convergence theorem in fixed bounded intervals (0, T), and we only need to check 
that the sequence {uh} is tight and uniformly integrable in time and space. Tightness 






in time is obvious in bounded intervals, and tightness in space is easily obtained 1 ^ 
Young’s inequality 


Vt > 0 : 


r?|xP^dx = 


uj^ ^ [uh\x\‘^Y dx < C 


gi(l/e)' 6Ll/e 


uniformly in h, t for some suitably small e > 0. Here we used the uniform bounds on 
the second moment m 2 {uh) < C and q — e < q < tq to control 

by local uniform bounds ||ith(^) ||L’’o(Rd) ^ Ct (obtained by propagation of initial 
integrability, Proposition 14.311 . The same propagation of integrability gives uniform 
bounds ||^t/i||Loo(o,T);L’’o(R'*) ^ C, thus by immediate interpolation we obtain 

equi-integrability in the form 

||^t/i||LP+e(0,r;L'2+£(R'*)) ^ 

for e > 0 suitably small (essentially such that 1 < g + e < tq). Applying Vitali’s con¬ 
vergence theorem gives strong convergence u/j —)■ u in L^([ 0 , T); L'^(]R'’*)) as desired 
and the proof is complete. 

□ 


We can now prove our main result. The proof of Theorem [3] is identical to that 
of Theorem [2] so we only establish the latter. 

Proof of Theorem\^ Step 1: convergence. Recall that Ufi{t),Vfi{t) G K, for all 
t,h, and that K, is L^-weak relatively compact. Using the approximate 1/2 Holder 
equicontinuity (15.7p and applying the previous refined Arzela-Ascoli theorem, we can 
extract a subsequence such that 

Vt > 0 : Uh{t),Vh{t) ^ u{t),v{t) in 

for some u,v £ C(0,T;P(M'^)) and u{t),v{t) G JC for all times. This entails the 
L°°(0, T; L”*(M'’*) nL^(M'’*)((l + |x|^)da:)) bounds. By standard truncation arguments 
we also get that u{t),v{t) are probability measures for all times, and because is 
l.s.c for the L^- weak convergence we can moreover take the limit in (|4.2I) to deduce 
that 1 1 —>■ u{t),v{t) are 1 / 2 -Holder continuous in {V,dw)- Since tt/i( 0 ), u/i( 0 ) = u^,v^ 
we can take the limit tt( 0 ),u( 0 ) = u^,v^, which together with u,v £ C^/^([ 0 , oo); P) 
shows that the limit u, v satisfies the initial condition at least in the sense of measures 
as desired. 

We claim now that 

^ in L^{6 ,T-LYrY) (5.9) 

for all 0 < (5 < r, and also 

G lY0,T-,lYrY)- 

To see this fix a test function ip £ L‘^{5,T', L°° {W^))y write for small h > 0 
f f {'Vu]f)ipdxdt = f f dxdt. 

Js JR<i Js JR<i 

By Corollary 14.11 and pointwise a.e convergence Uh ^ u we can assume ^ 

in L^((5, T; L^(M'’*)) for fixed 0 < 5 < T, and by diagonal extraction we can 
assume that the limit is independent of 5,T. By (15.2p with (7 = m < tq it is 



easy to get in L^((5, T; L^(M'^)). As a consequence we can pass to 

the limit 



{Vu^)pdxdt 2 
n\0 



dxdt 



Vu^pdxdt 


to obtain (15.9p . In particular by Corollary 14.II we see that VO < (5 < T, it holds 


Vu™ 




L2(5,r;Ll(K‘')) 


< 2 lim inf 11 
h\o " 


L°°(5,T;L2(Krf)) 


Vu 


m/2 

h 


L2(5,T;L2(Md)) 


< C (1 + T )^/2 


uniformly in 5 > 0, whence Vu^ G L^(0, T; for all T > 0. 

For the drift terms, recall from m that VC/(x),VC(x) are at least locally 
bounded. From Theorem [5] we have Uh,Vh —t u,v at least in L^(0, T; 
thus Uh^U,Vh^y uVU,v\/V in L^(0, oo; when tested with compactly 

supported functions p{x). Moreover by uniform bounds on the second moments 
and linear behavior of VC/, W it is easy to check that the limit uS/U, v'VV G 
L°°{0,oo-,L^(R^)). 

Regarding now the coupling terms rt/jV'i/’ft, u/iV'0/i, note from (|5.2p with q = 
2dl (d + 1) < xq that we have in particular Uh, Vh ^ u,v in LP{6, T ; fQj- 

all p G [l,oo). By strong j^ 2 d/(d-i) continuity in pHLS-.Sp we 

thus obtain V'lph = {VG)*[uh-Vh] ^ (VG)*[u-u] in LP(d,T; 

for all p G [1, oo), so by Holder inequality 


Uft VhV^Ph ^ vV^P in 1^(6, T; (M'")) 


for all p G [1, oo) and 0 < <5 < T. Using the L°°{5, T] bounds for u^, Vh 

this gives LP((5, T;L^(M'^)) bounds uniformly in p > 1 and <5 > 0, thus uV'tp,v'V'ip G 
L°°(0,T;Li(R'=')) C L^{0,T-L^{R<^)) for all T > 0. 

Step 2: the weak solution. Fix any test-function <p G and 0 < Ti < T 2 , 

and let A^i = [Ti/h\,N 2 = [^ 2 /^]. Let be the optimal map in = 

{y^u/‘\ Expanding 


p{x) - p(Vg(")(x)) = (V<p(Vg(")(x))) • (x - Vq^^\x)) 

+ 0(\x - Vg(-)(x)|2||Z/VllLoo(K.)). (5.10) 




Taking (^ = Vy? in the Euler-Lagrange equation (|4.5I) . and summing from n = Ni 
n = N 2 — 1 we compute by (|4.1I) . (|4.5p and (15.101) that 


[ {uh{T 2 )-Uh{Ti))ifdx 
Jr<^ 




)(/?dx 


N2-1 


E /,(• 


(n+l) (n) 


- Ul 


)(/J dx 


n=Ni 

N2-I 


/ \ipo — ip\ dx 

»=iVi 

V2-1 ^ 

- Id, V(/? O Vjr^ dx 


n=A^i ' 


iV2-l 


+ 




n=A^i 


A^2-1 


= yy h / Vv:>) - 4"+'^ 

n=iVi ^ 

Integrating by parts and exploiting (|4.ip . at the continuous level this becomes 


h 


dx 


C’(/i||^VllL<-(Kd)) + [ {uh(T2)-Uh{Ti))ipdx 

JR^ 

pN2h p 

= -/ / (VK(t))™,V<^) + (V17,V^K(t) + u;,(t)(VV’/.(t),V<^) 

JNih JRd- L 


f'N2h 
iNih 

By step 1 we can take h —?> 0 as 


dxdt. 



— u(Ti))</?dx 



, V(/j) + u{VU, V(/?) + uCVip, V(/j)) dxdt 


to obtain (12.8p . The equation for v is similarly obtained. 

Step 3: Energy bounds and further regularity. We first establish the energy 
bound (|2.1ip . Arguing as in step 1, using p5.2p with q = 2d/{d+2) < tq, and (IHLS-3P 
it is easy to conclude that S/iph —?> S/ip in 00 ; for all p € [l,oo). In 

particular Vd>h{t) —?> V'I'(t) in and 


£cpi{uh{t),Vh{t)) £cp\{u{t),v{t)) a.e. t > 0. 

From the energy control ||Vt/’/i|| L°°( 0 ,oo;l, 2 (]Rd)) < C we also have V-^ G L°°(0, 00 ; L^(M'^)). 
Similarly, from the 00 ; with g = m < rg it is easy to get 

£6is{uh{t),Vh{t)) £Ais{u{t),v{t)) a.e. t > 0. 












I^or the potential energy, we have 0 < Uh{t,x)U{x) u{t,x)U{x) a.e. x G for 
a.e. fixed t > 0, and similarly for VhV- By Patou’s lemma we conclude that 

£e^tiu{t),v{t)) < \immf£ej^t{uh{t),Vhit)) a.e. t > 0. 
h^O 

Summing S = f’difr + Sext + ^cpi with £{uh{t),Vh{t)) < £{vP,v^) at the discrete 
level finally entails the desired energy bound (12.lip . Note that the potential energy 
fext is the most problematic term, because we cannot a priori conclude equality 
£extiu{t),v{t)) = lim£ext{uh{'t):Vh{t)) in the last display (in fact strong convergence 
would imply convergence of the second moments, thus convergence in dw)- As a 
consequence we only retrieve the one-sided inequality, and we are unable to conclude 
that the total energy is monotone nonincreasing in the limit. 

Turning now to the propagation of initial regularity (I2.12p , assume that the initial 
datum G for some p G [l,oo]. If p < oo then Proposition 14.31 bounds 

Uh{t),Vh{t) in uniformly in h with exponential control 

SUp^ ^ + |b°||LP(Ed)^ . 

Up to extraction of a further subsequence we can assume that Uh,Vh ^ u,v in 
L“^([0, oo); LP{W^)), which immediately gives (12.121) . If now u^, G clearly 

(12.121) holds for arbitrarily large p. Our claim then easily follows by letting p ^ oo 
and the proof is achieved. 

□ 


6 . Appendix 

For p > 1 let be the weak-L^’ spaces, which coincide with the usual Lorentz 

space The natural Banach norm is 

lkllLS,(R0 = (6.1) 

where w*{t) is the symmetric-decreasing rearrangement of w{x). 

Proposition 6.1. Denoting $ = (—A)“^r(; = G *w, the Dirichlet energy 

w G L^(M'’*) !->• £d{w) = J |V<l>p dx G [0, -|-oo] 
is lower semi-continuous for weak convergence. 

Proof. Let Wn ^ rc in If liminf ££i{wn) = +cx3 our statement is trivial, so up 

to extraction of a subsequence we may assume that liminf £’£)(tc„) = \\va.£D{wn) = 
C < -|-oo, in particular we have that 

lim ||V<l'„||^ 2 (Rd) = liminf £(£)(u;„) < -|-oo. (6-2) 

Now since Wn w in we see that Wn is bounded in L^(M'’*), hence bv (IHLS-2l) . 

||^n||^d/(d-2)^^d) < C'lknllLl(Rd) < C. 

Since = (L'’*/^’^(M'’*))^ is a topological dual we can 

also assume, by the Banach-Alaoglu theorem, and up to a further subsequence, that 

in 

By (|6.2I) and up to a subsequence we see that 

||V^'||^2(Rd) < liminf IIV$n|li2(Rd) = liminf ^(^(rcn). 












As a consequence it suffices to prove that = G*w, since then ££i{w) = 
lim inf £D{wn)- 

Set ^ = G * w ^ and let us prove that <h — $ = 0. Since — A<f>„ = 

Wn ^ w = — A^> in we have in particular that $ — $ is harmonic. Because 

harmonic tempered distributions are polynomials and C we get 

that (M'^) is polynomial. By (|6.ip we see that the polynomial $ — $ 

decays at infinity, hence <f> — $ = 0 as claimed and the proof is complete. □ 
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